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' Abstract 



This paper is addressed to a study of the nuh controhabihty for the semihnear parabohc 
] equation with a complex principal part. For this purpose, we establish a key weighted 

identity for partial differential operators (a + i(3)dt + J2 (with real func- 

j,k=i 

' tions a and by which we develop a universal approach, based on global Carle- 

. man estimate, to deduce not only the desired explicit observability estimate for the 

\ linearized complex Ginzburg-Landau equation, but also all the known controllabil- 

ity/observability results for the parabolic, hyperbolic, Schrodinger and plate equations 
that are derived via Carleman estimates. 
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95 ; 1 Introduction and main results 

> . 
•1— I . 

^ . Given T > and a bounded domain Vt of IR" (n G N) with C boundary F. Fix an open 
^ ! non-empty subset a; of f2 and denote by Xlo the characteristic function of u). Let ojq be 
another non-empty open subset of such that uJq d u. Put 

Q=(0,r)xfi, S = (0,T)xF, Qo = (0,T) X 6^0- 

In the sequel, we will use the notation yj = i/x , where Xj is the j-th coordinate of a generic 
point X = {xi, ■ ■ ■ , Xn) in R". In a similar manner, we use the notation zj, vj, etc. for 
the partial derivatives of z and v with respect to Xj. Throughout this paper, we will use 
C = C{T, to denote generic positive constants which may vary from line to line (unless 
otherwise stated). For any c G(D, we denote its complex conjugate by c. 
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Fix a^'=(-) e Ci'2(g; R) satisfying 

a^f'^t, x) = a''\t, x), {t, x) eQ, j,k^l, 2, 
and for some constant Sq > 0, 



Next, we fix a function /(•) e C^((D) satisfying /(O) = and tlie following condition: 

1/(^)1 



lim 

s— >oo 



lnV2 



0. 



(1.1) 
(1.2) 

(1.3) 



Note that /(•) in the above can have a superlinear growth. We are interested in the following 
semilinear parabolic equation with a complex principal part: 



;i + tb)yt - {o.^^yj)k = XujU + f{y) in Q, 
j,k=i 



[ 1/(0) = yo 



on E, 
in Q, 



(1.4) 



where i = & G IR. In (1.4), y = y{t, x) is the state and u = u{t, x) is the control. 

One of our main objects in this paper is to study the null controllability of system (1.4), 
by which we mean that, for any given initial state yo, find (if possible) a control u such that 
the weak solution of (1.4) satisfies y{T) — 0. 

In order to derive the null controllability of (1.4), by means of the well-known duality 
argument (see [16, p. 282, Lemma 2.4], for example), one needs to consider the following 
dual system of the linearized system of (1.4) (which can be regarded as a linearized complex 
Ginzburg-Landau equation): 

Gz — q{t, x)z in Q, 

z^O on E, (1.5) 

z{T) = zt in Q, 
where g(-) G L°°(0, T; L"(fi)) is a potential and 



j,k=i 



(1.6) 



By means of global Carleman inequality, we shall establish the following explicit observ- 
ability estimate for solutions of system (1.5). 

Theorem 1.1 Let a^^{-) e C^'2(Q; R) satisfy (1.1)-(1.2), q{-) e L°°(0, T; L"(^])) and b e 
R. Then there is a constant C > such that for all solutions of system (1.5), it holds 



k(0)|L2(n) < C(r)|z|i2((o,T) 



where 



C{r) = Coe^""^', Co = C(l + b'), r = |g|Loc(o,T;L«(Q)). 
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(1.7) 
(1.8) 



Thanks to the dual argument and the fixed point technique, Theorem 1.1 imphes the 
following controllability result for system (1.4). 



Theorem 1.2 Let a^^{-) e C^'^iQ; H) satisfy (1.1)^(1.2), /(•) e Ci((D) satisfy f{0) = 

and (1.3), and 6 G R. Then for any given yo G L'^{^), there is a control u G -^>^((0, T) x uj) 
such that the weak solution y{-) G C([0, T]; L^(0)) fl -^^(0, T; ifo(^)) of system (1.4) satisfies 
y{T) = zn n. 

The controllability problem for system (1.4) with 6=0 (i.e., linear and semilinear 
parabolic equations) has been studied by many authors and it is now well-understood. 
Among them, let us mention [4, 5, 13, 9] on what concerns null controllability, [7, 8, 9, 28] for 
approximate controllability, and especially [29] for recent survey in this respect. However, 
for the case 6 7^ 0, very little is know in the previous literature. To the best of our knowledge, 
[10] is the only paper addressing the global controllability for multidimensional system (1.4). 
We refer to [20] for a recent interesting result on local controllability of semilinear complex 
Ginzburg-Landau equation. 

We remark that condition (1.3) is not sharp. Indeed, following [4], one can establish the 
null controllability of system (1.4) when the nonlinearity f{y) is replaced by a more general 
form of f{y, Vy) under the assumptions that /(•, •) G C^((D^+"), /(0, 0) = and 



where p = (pi, • • • ,Pn)- Moreover, following [5], one can show that the assumptions on the 
growth of the nonlinearity /(y, Vy) in (1.9) are sharp in some sense. Since the techniques 
are very similar to [4, 5], we shall not give the details here. 

Instead, as a byproduct of the fundamental identity estabhshed in this work (to show 
the observability inequality (1-7)), we shall develop a universal approach for controllabil- 
ity/observability problems governed by partial differential equations (PDEs for short), which 
is the second main object of this paper. The study of controllability/observability problem 
for PDEs began in the 1960s, for which various techniques have been developed in the 
last decades ([1, 3, 13, 17, 29]). It is well-known that the controllability of PDEs depends 
strongly on the nature of the system, say time reversibility or not. Typical examples are 
the wave and heat equations. It is clear now that there exists essential differences between 
the controllability/observability theories for these two equations. Naturally, one expects to 
know whether there are some relationship between these two systems of different nature. 
Especially, it would be quite interesting to estabhsh a unified controUabihty /observability 
theory for parabolic and hyperbolic equations. This problem was posed by D.L. Russell in 
[21], where one can also find some preliminary result; further results are available in [18, 26]. 
In [15], the authors analyzed the controllability/observability problems for PDEs from the 
point of view of methodology. It is well known that these problems may be reduced to 
the obtention of suitable observability inequalities for the underlying homogeneous systems. 








|(s,p)Hoo ln^/2(l + |s| + \p\) 



(1.9) 



lim 

|(s,p)|^oo 
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However, the techniques that have been developed to obtain such estimates depend heavily 
on the nature of the equations. In the context of the wave equation one may use multi- 
pliers ([17]) or microfocal analysis ([1]); while, in the context of heat equations, one uses 
Carleman estimates ([8, 13]). Carleman estimates can also be used to obtain observabihty 
inequalities for the wave equation ([25]). However, the Carleman estimate that has been 
developed up to now to establish observability inequalities of PDEs depend heavily on the 
nature of the equations, and therefore a unified Carleman estimate for those two equations 
has not been developed before. In this paper, we present a point-wise weighted identity for 

partial differential operators (a -|- iP)dt + J2 dk{a^^dj) (with real functions a and /3), by 

j,k=i 

which we develop a unified approach, based on global Carleman estimate, to deduce not 
only the controllability/observability results for systems (1.4) and (1.5), but also all the 
known controllability /observability results for the parabolic, hyperbolic, Schrodinger and 
plate equations that are derived via Carleman estimates (see Section 2 for more details). We 
point out that this identity has other interesting applications, say, in [19] it is applied to 
derive an asymptotic formula of reconstructing the initial state for a Kirchhoff plate equation 
with a logarithmical convergence rate for smooth data; while in [11] it is applied to estab- 
lish sharp logarithmic decay rate for general hyperbolic equations with damping localized in 
arbitrarily small set by means of an approach which is different from that in [2]. 

The rest of this paper is organized as follows. In Section 2, we establish a fundamental 
point-wise weighted identity for partial differential operators of second order and give some 
of its applications. In Section 3, we derive a modified point-wise inequality for the parabolic 
operator. This estimate will play a key role when we estabhsh in Section 4 a global Carleman 
estimate for the parabolic equation with a complex principal part. Finally, we will prove our 
main results in Section 5. 

2 A weighted identity for partial differential operators 
and its applications 

In this section, we will establish a point-wise weighted identity for partial differential opera- 
tors of second order with a complex principal part, which has an independent interest. First, 
we introduce the following second order operator: 

Vz^{a + iP)zt + Yl {a'''zj)k, m e IN. (2.1) 
j,k=i 

We have the following fundamental identity. 

Theorem 2.1 Let a, (3 e C2(Ri+™;R) and a^'^ G Ci'2(Ri+'"; R) satisfy a^^ = a^^ {j,k = 
1,2, •••,m). Let z, v e C^R^^""; (E) and £ e C2(Ri+"^;R). Set 6 = and v = 9z. 
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Then 

e(Vzh + Vzh) + Mt + div V 

m 1 



j,k,j',k'=l 



+i J2 [{Pa^'"ij)t + a^''{Pit)j]{vkV - Vkv) - ^ a^^ak{vjVt + VjVt) 
j,k=i j,k=i 



+i 
where 



j,k=i 



j,k=l j,k=l 



and 



M 



A 



j,k=i 



j,k=i 



iV — VjV j 



— aa'^{vjVt + VjVt) 



(2.2) 



(2.3) 



(2.4) 



V ^[V\---,V\---,V'^l 

m 

3,3',k'=l 

^(^a'^'a^'^ - a^''a^'''')£j(vfVk' + Vj>Vk') - '^a^^{vjV + Vjv) 

A ^ 

B = {aHt)t + (/3'4)t - (aA)* - 2[ ^ + 

m m 
j,k=l j,k=l 

Several remarks are in order. 

Remark 2.1 We see that only the symmetry condition of {a^^)mxm is assumed in the above. 
Therefore, Theorem 2.1 is applicable to ultra-hyperbolic or ultra-parabolic differential opera- 
tors. 

m 

Remark 2.2 Note that when = — E i^^^^j)k (^nd a{t, x) = a, P{t, x) = b with a, 6 e 

j,k=i 

Theorem 2.1 is reduced to [10, Theorem 1.1]. Here, we add an auxiliary function \l/ in the 
right-hand side of the multiplier Ii so that the corresponding identity coincides with [12, 
Theorem, 4-1] for the case of hyperbolic operators. Moreover, we will see that the modified 
identity is better than [10] in some sense. 
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Remcirk 2.3 By choosing a{t,x) = 1, l3{t,x) = and^ = ~2 ^ {a^''£j)k in Theorem 2.1, 

j,k=i 

one obtains a weighted identity for the parabolic operator. By this and following [22], one 
may recover all the controllability/observability results for the parabolic equations in [4] and 
[13]. 

Remark 2.4 By choosing aP^(t,x) = a^^{x) and a(t,x) = (3{t,x) = in Theorem 2.1, 
one obtains the key identity derived in [12] for the controllability/observability results on the 
general hyperbolic equations. 

Remark 2.5 By choosing {a^^)i<j^k<m to be the identity matrix, a{t,x) = 0, (3{t,x) = 1 
and ^' = —Ai in Theorem 2.1, one obtains the pointwise identity derived in [14] for the 
observability results for the nonconservative Schrddinger equations. Also, this yields the 
controllability/observability results in [27] for the plate equations and the results for inverse 
problem for the Schrddinger equation in [6]. 

Remcirk 2.6 By choosing {oP^)i<j,k<m to be the identity matrix, a{t.,x) = 0, [3{t,x) = 
p{x) and — —Ai in Theorem 2.1, one obtains the pointwise identity for the Schrddinger 
operator: ip{x)dt + A. Further, by choosing 

i{t, x) =sif, ip = eT(l^-^ol'-^l*-*ol') 

with 7 > 0, c > 0, .Xq G \ Q and assuming Vlogp ■ {x — xq) > —2. One can recover the 
fundamental Carleman estimate for Schrddinger operator ip{x)dt + A derived in [24, Lemma 
2.1]. 

Proof of Theorem 2.1. The proof is divided into several steps. 

Step 1. By 6 = 6^, v = 9z, we have (recaUing (2.4) for the definition of /i) 

m 

eVz = (a + i(5)vt - (a + i(5)itv + 

3,k=l 



j,k=l j,k=l j,k=l 



(2.5) 



where 



/2 = avt - i(3itv d^^^jVk + '^v. (2.6) 

3,k=l 

Hence, by recaUing (2.3) for the definition of /i, we have 

e{Vzh + V^h) = 2|/i|' + {IJ2 + hli). (2.7) 

Step 2. Let us compute hi 2 + hh- Denote the four terms in the right hand side of Ii 
and I2 by If and respectively, j — 1, 2, 3,4. Then 

+ + II) = -{(3Ht\vW + {(3Ht)t\v\\ (2.8) 
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Note that 



Hence, we get 



2vvt = {\v\'^)t - {vvt - vvt), 
2vvk = - {vvk - vvk). 



mil + li) + mil + II) 

m 

= -2i 53 [{Pa^%vvk)t - {pa^'ij)tvvk 



j,k=i 



rn 

+2i ^ [{(3a^%vvt)k - {Pa^''ij)kvvt] - i(3^{vvt - Vtv) 



j,k=i 



j,k=i 

m 



j,k=l 



Next, 



j,k=l j,k=l 

mt +t+t+ It) + till + ii + ii + It) 

m 

j,k=i 

m 

+ {c^Ht)t\v\'' -2[Y {aa^'Wk + «^A]br- 
j,k=i 

Noting that a^'^ = a''\ we have 

mn+ii)+mn+ii) 

m mm 



j,k=i 

m 



j,k=l 



j,k=l 



j,k=l j,k=l 



VjV - VjV) 



+ t Y <i"'W^t)ki 



VjV — V- 



j,k=l 



rn 

Y [aa^'^ivjVt + VjVt) + if5a^''it{vjV - Vj 



j,k=i 

m m 

- J2 ("a^S'^fc)*- Yl 0,^''0Ck{VjVt + VjVt) 
j,k=l j,k=l 
2 m m 

+ 2 J2 (^^^%i'"j'^k + VkVj) + i Y ci^^{(3Qk{vjV - vyv). 

j,k=l j,k=l 



Using the symmetry condition of a^^ again, we obtain 



j,k,j',k'=l 

in 

j,k,j',k'=l 



(2.12) 



j,k,j',k'=l 



By(2.12), we get 



If /| + If I? 



rib in 



j,k,j',k'=\ 

m 



j,k,j',k'=l 



(2.13) 



j,k,j',k'=l 



j,k,j',k'=l 



Further, 



III fit, 



j,k=i 

m 



j,k=l 



j,k=i 



(2.14) 



Finally, 



j,k=i 

+ t + t + 1) + + 4^ + II + Z^^) 
= {aA\vW - {aA\\v\'' - 2 E (^''^^kH^ 



+2[ 5:(a^%A), + A* 



\v?. 



(2.15) 



.Step 5. By (2.8)-(2.15), combining all '^-terms', all '^-terms', and (2.7), we arrive at 
the desired identity (2.2). □ 

We have the following pointwise estimate for the complex parabolic operator Gz. 

Corollary 2.1 LetheB. and a^''{t,x) e C^'2(Il^+"; R) satisfy condition (1.1). Let z, v E 
C^iB}-^"^; (D) and ^, ie C2(]El^+"^; R). Set 9 ^ and v ^ Bz. Put 



* = -2 ^ {a^%),. 
j,k=i 



(2.16) 



Then 



e'\Gz\' + Mt + divV 



n 1 



j,k,j',k'=l 



(2.17) 



+i6 ^ (af + 2a^%t)(l;fe^; - ^;feTJ) - i6 E («''^j)ik(^^t - ^^t) + ^l^^l", 

j,k=l j,k=l 



where 



M 



a>^{vjVt + VjVt) 



[VjV-VjV) 



(2.18) 



j,f,k'=l 

+ {2a'^' a>''^ - tt>^a^'^'^lj{vj,Vk' +Vj'Vk') - '^a^^{vjV + Vjv) 
+a^''{2A£j + - 2£/t)|w|2}, 

n 

B={1 + e)£u - A - 2 [ 5: {a^H,£t)k + 

j,k=i 

n n 

Proof. Recalling (1.6) for the definition of Gz, taking m = n, a{x) = 1, P{x) = 6 in 
Theorem 2.1, by using Holder inequality and simple computation, we immediately obtain 
the desired result. □ 



3 A modified point-wise estimate 

n 

Note that the term ib ^ {a^''ij)k{vvt — vvt) in the right-hand side of (2.17) is not good. 

j,k=i 

In this section, we make some modification on this term and derive the following modified 
point-wise inequality for the parabolic operator with a complex principal. 

Theorem 3.1 Let h e and a^^{t,x) E Ci'2(]Ri+"; R) satisfy (1.1)-(1.2). Let z, v e 
C2(1R^+"^; (D) and ^ G C2(1R^+"^; R) satisfy (2.16). Put 



e^e^, 9z. (3.1) 



Then 



2e^\Gz\'^ + Mt + diyV 

n 

> E c'\vkVj + VkVj) + B\ 



(3.2) 

-1 n 

7,fc=l + i',fe' = l 



j,k=l 

n 
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where M, V^, B is given by (2.18) and 

ib 



j,j',k'=l j,k,j',k'=l 



1 + 62 



1 + 62 



d'^\z 



n 11 

ci'= E [2{a^'%'ha^''-ai';a^'''i, + -ai' + ^^^^ 



j',k'=i 
B^B 



262 
b^ 



j,k,j',k'=l j,k=l 



Y: {2a^%{a^'''ir)k'j + {a^'{a^'''if)k'j)A. 

j,k,j',k'=l 

Proof. We divided the proof into several steps. 
Step 1. Note that v — 9z, it is easy to check that 

VVt - VVt = 9'^izZt - ZZt), VVk - VVk = 0'^{zZk - ZZk). 

RecaUing (1.6) for the definition of Gz and by (3.4), we have 

n n 

-ib ia^''ij)k{vvt - VVt) = -ib {a^''^j)kO^{zzt - zzt) 

j,k=l j,k=l 

ji.n2 n 

Y (1 - ib)Gzz - (1 - ib)Gzz 



1 + 

— - Y E [ia'%)kZ-{a^''zj)kz) 

j',k'=l j,k=l 



2n2 n 



E i^''''^f)k' E li<^'%)kz + {a^'z,: 
j',k'=i j,k=i 

3 



kZ, 



k=l 



Step 2. Let us estimate Jk {k = 1,2,3) respectively. 
First, note that v — 9z, we have 



Ji = 



> 



•L/i2 n 

Y {a^'^Qk [(1 - ib)Gzz - (1 - ib)G 



1 + 62 



zz 



j,k=l 



(1 - ib)9Gz 


2 


VI + 62 





= -9^\Gz\' 



n 2 
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Next, by using (3.4) again, we have 



J2 = 



i',fc'=i i,fc=i 



ib 



1 + 62 



^ [6*2 {a^''''ij,)k'a^''{zjz - Zjz) 



j,k,j',k'=l 



lb 



E [24(o^'''^,')ik' + J a^'izjZ - Zjz) 



1 + 62 

j.k.j' ,k'=l 



j,k,j',k'=l 

ib 



VjV) 



+ 



lb 



1 + 62 



j,k,j',k'=l 

Next, noting that a^'^ satisfy (1.1), and recaUing v — 9z, it follows 

62^2 n 



T-TJ E ^r)k' E [(a^%).^+(a^'=^,).^) 
52 " 

1 + ^,2 E { 

"-^^ j,k,j',k'=l 



^^{a^''''ij')k'a^''{zjZ + Zjz) 



{a^'^' ^j')k'a^^{zjZk + ZjZk) 



1 + 62 



j,k,j',k'=l 



E \\0'{a^''''ir)k'] \z\ 



eHa^''''ii')k'a^''iziZ + 





r+62 ^ "i- 

962 n 

+7-7^ E [(a^'''^/).'a^''4| 



^'(a^'''4)fe' 



j,k,j 

62 



1 + 62 



E (a^''''^j')fc'a^''(^Jj^;ifc + -yjt^ifc) + {2 E a^'^^Aia^''"' ij')k 

j,k,j',k'=l j,k,j',k'=l 
n n 

+2 a^%{a''''ej')k'j+ E («''(«''''^i')A.'.)Jl^r 

j,k,j',k'=l j,k,j',k'=l 
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where we have used the following fact: 

n 

j,k,j',k'=l 

n n 

= {a^''''£r)k'a^''{vjVk + VjVk) + 2 ^ {a^'''' ej,)k'a^''ijek\v\ 

j,k,j' ,k'=l j,k,j' ,k'=l 

n n 



j,k,j',k'=l 



j,k,j',k'=l 



and 



E {[9\a^''''i,),,].a^'}jz\' 

j,k,j',k'=l 
n 

= Y {e^[2a^%{a^''''ir)k' + a^''{a^''''ij')k'j]}jz\^ 

j,k,j',k'=l 

n n 

j,k,j',k'=l j,k,j',k'=l 



+2| E 



+ E {a^\a^'''er)k'j)^}\^ 

j,k,j',k'=l 



Step 3. Noting that a^^ satisfying (1.1)- (1-2), we have the following fact. 
1 " ■ 6 " ■ 



1 



1 + 62 



E 0'''^ ^jVk — 'iibikivjV — Vjv) + Ab'^ij£k\v\ 



j,k=i 



Ab 



2 n 



1 + 62 



E a'%4bl 



1 ™ n 

E a^*=(T;, - 2z6^,^;)(^;fc - 2^64^^) - jip^ E 



j,k=i 

,,|2 



1 + 62 



> 



462 " 



1 + 62 



E «'%4|^^|'. 



(3.9) 



(3.10) 



(3.11) 



Finally, combining (2.17), (3.4)-(3.8) and (3.11), we arrive at the desired result (3.2). 



4 Global Carleman estimate for parabolic operators 
with complex principal part 

We begin with the following known result. 

Lemma 4.1 ([13], [23]) There is a real function ip £ C'^{Vt) such that ip > infl and ip = 
on dQ and \Vtp{x) \ > for all x E Q \ uq. 
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For any (large) parameters A > 1 and fJ> > 1, put 

i^Xp, <p{t,x)^j^^, p{t,x)^ . (4.1) 

For /c = 1, • • • , n, it is easy to check that 

£t = Xpt, ij = X/j^ifiipj, Ijk = XjJi^ipiijipk + >^IJ-vi'jk (4.2) 

and 

In the sequel, for A; e IN, we denote by 0{p^) a function of order ji^ for large ii (which is 
independent of A); by 0^j,{\^) a function of order A'^ for fixed ii and for large A. 

We have the following Carleman estimate for the differential operator G defined in (1.6): 

Theorem 4.1 Let 6 e R and a^^ satisfy (1.1)-(1.2). Then there is a po > such that for 
all p > Pq, one can find two constants C — C{p) > and Ai = Xi{p), such that for all 
z e C([0,T];L2(Q))nC((0,r];//o^(fi)) and for allX>Xi, it holds 



A V / (p^O^\z\^dtdx + Xp^ [ <f9^\Vz\^dtdx 

< C{1 + b^) \ f e^\Gz\^dtdx + AV / '/d^\z\^dtdx 

'-Jq J{0,T)xui 

Proof. The proof is long, we divided it into several steps. 

Step 1. Recalling (3.3) for the definition of c^'', by (4.2)-(4.3), we have 



d {vkVj + VkVj) 

j,k=i 

n 

j,k,j',k'=l 



(4.4) 



1 1 



2 * 2(1 + 62) ^ ^ 



{vkVj + VkVj) (4.5) 



j,k=i 

+^7^— T^W a^''a^''''iJj'ipk'{vkVj + VkVj) + X(pO{p)\Vv\'^. 

y ' j,k,j',k'=l 

On the other hand, by (4.2)-(4.3), recalling (2.16) and (2.3) for the definitions of ^ and A, 
it is easy to check that 

n n 

* = -2A/iV E ^''Mk + X^O{p), A = A^/zV' E (''"^j^k + V^'O^(A). (4.6) 

j,k=l j,k=l 
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Next, recalling (2.3) and (2.18) for the definition of A and B, respectively. By (4.2)-(4.3) 
and combining (4.6), we have 



7t IL 

B ^2J2 (^''^jAk -2AY: 

j,k=l j,k=l 

TL fL Tl 

+ (1 + b')iu - A - 2 X: a^%etk + 2£t E {a^'ij)k + E 

j,k=l j,k=l j,k=l 

= 2aW| E +aVc(/x') + ¥''c)^(a2). 

j,k=i 

Hence, by recaUing (3.3) for the definition of B, we have 

~ 2 I " 2 

j,k=l 



1 + 62 

^. By (4.2)-(4.3), we have 

n 

\ib E («f + 2a^%t)(Ufe^; - Vkv) 
j,k=i 



<C\^ip^ E (^^^i^ji^kV — Vkv) < CA/Lt(^| E oP^'i'jVk +CXfj,(p 

j,k=l j,k=l 

It is easy to see that (4.9) can be absorbed by (4.5) and (4.8). 
Similarly, by using (4.2)-(4.3) again, we have 

ib 



— — E {a^''''£j')k'ja'\vkV - Vkv) 
j,k,j',k'=i 



< 



1 



< 



1 + 62 

1 



1 + 62 

C 

+ 



j',k'=l 



j',k'=i 



< 



1 + 62 

1 



bXi/ip^ E a^^i^jivkV - Vkv) 

j,k=i 



1 + 62 



A/iV E ^^''i'j'^k 



j,k=i 



1 + 62 



A/iV-1 E a^'V'iV'fe 



+ 



c 



1 + 62 



C62 



1 + 62' 



Therefore (4.10) also can be absorbed by (4.5) and (4.8). 



(4.7) 



(4.8) 



(4.9) 



(4.10) 
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Combining (4.5), (4.8)-(4.10), by (1.2), we end up with 



The right-hand side of (3.2) 
1 



> 



1 + 62 

2 



0',k'=l j,k=l 



+ 



> 



1 + 62 

1 



1 + 62 

2 



j,k=i 



(4.11) 



+ 



1 + 62 



3,^3/ 



Step 3. Integrating (3.2) on Q, by (4.11), noting that ^(0) = e{T) = 0, we have 



A/^='|V^|2 + A0(/x) 



Wv\ dxdt 



+ / [aVI VzM^ + X^Oiijc") + 0^(A2)1 



(4.12) 



< C 



\eGz\ 



LHQ) 



+ / div V ■ vdxdt 
JQ 



Next, recalhng (2.4) and (3.3) for the definition of V and V . Noting that = and 

dv J 
Si/' 



Vi = (which follows from (1.5) and ^Is = 0, respectively), by (4.2) and Lemma 4.1, we 



have 



/ div V • udxdt — / div V • vdxdt 
JQ JQ 



V^V^I^ll^f ( E a^^i^kfdtdx < 0. 



dv |2 



(4.13) 



j,k=i 



On the other hand. 
Left-side of (4.12) 



> 



ip(Xi^^\Vij\^ + XO{i^))\Vv 
+^^(aV|VV'|' + X'Oifi') + ^'0,iX'))\v\ 

+(^=^(aV|vV'|' + x'oii,') + ^'o,ix') 

-CX^i^ [ (p{\Vv\^ + X^iJ,\''\v\^)dtdx. 

jQn 



dtdx 



dtdx 



(4.14) 
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By the choose of we know that min jV^'l > 0. Hence, there is a //q > such that 
for all fJ,> /Iq, one can find a constant Ai = Ai(/u) so that for any A > Ai, it holds 



r [ ip\XtJ,^\Wij\^ + XO{fj,)]\Vv\^dtdx 
Jo Jn\ujo 

+ / ^3[AV|V^r + A30(;.3) + ^3o^(A2) 

> coXfi^ r [ ^{\Vvf + X'^^^ip'^\v\'^)dtdx, 
Jo Jn\u)o 



\v\ dtdx 



iCl\ujo 

where cq — min ( min |V^|^, min |VV'|^) is a positive constant. 
Next, note that v = 9z, by (4.2), we have 

Zj — 0~^{vj — £jv) — 0~^{vj — Xfj,(pipjv), Vj — 9{zj + ijz) — 9{zj + Xfxipijjjz). 
By (4.16), we get 

^eWvz\^ + AVVl'^l') < I v^l' + AW>r < C9\\vz\^ + aVV'I-^I')- 

Therefore, it follows from (4.14)-(4.15) and (4.17) that 

X// f ^9^{\Vzf + X'^i/(p^\z\^)dtdx 

JQ 

= A// r ( [ + [ )ipe\\Vz\^ + X^fi^ip'\z\^)dtdx 
Jo ^ Jn\ujo Jujo ^ 

dtdx 

J Q 

+ f (^'[AV|VV'|^ + A30(/x3) + (^30^(A2)]|^;|'dtda; 

+X^^ [ ip9WVz\'' + X''iJ,\''\z\^)dtdx}. 
JQo 

Now, combining (4.12)-(4.13) and (4.18), we end up with 

A//2 / ip9WVz\^ + X^iJ,\W)dtdx 
Jo 



< C 



9'\Gz\'dtdx + A/x' / <^9WVz\' + X' ii'if'\z\')dtdx 

JQo 



(4.15) 



(4.16) 



(4.17) 



(4.18) 



(4.19) 



Step 4- Finally, we choose a cut-off function ( G C^{u; [0, 1]) so that ^ = 1 on a;o. Then 



[C^9'\{l + ib); 



C\{1 + ib)z\'{<p9')t + CV^'(1 + ib){l - ih){zzt + zzt). (4.20) 
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dtdx 



By (1.6), (4.20) and noting ^(0, x) = e{T, x) = 0, we find 

= / C\\{l + ib)z\\ip9% + ipe^{l + ib){l-tb){zzt + zzt) 
jQo 

= / (:^e^\{l + ib)z\\cpt + 2Xcppt) + ^{l-ib)z\-J2{a"'z,)k + Gz 
+ip{l + ib)z[ - + G^] } 

= / e^{C\{i + ib)z\''{ipt + 2\ippt) 



+CVI^a^'''[(l - ih)zjZk + (1 + ib)zjZk 



X] [(1 ~ ib)zzjiik + (1 + ib)zzji)k 
+2Xii('^(p'^ a^^ [(1 - ib)zzjipk + (1 + ib)zzj'ijjk 
+2Cip Yl (^^^ [(1 - ib)zZjC,k + (1 + ih)zZjC,k 

+CV[(1 - ib)zGz + (1 + z6)zG^]}dic/x. 
Hence, by (4.21) and (1.2), we conclude that, for some 5 > 0, 
2 f CfO'^lVzl'^dxdt 

JQo 

= / eHC\{i + tb)z\\^t + 2\^pt) 

JQo 

+fiC'^(p Y [(1 ~ ib)zzjtpk + (1 + ib)zzjipk 
+2A//C^(/7^ Y 0'^'' [(1 - ib)zzj'ipk + (1 + ib)zzjipk 

j,k 

+2C^ Y a'" [(1 - ib)zzjCk + (1 + ib)zZjCk 

j,k 

+CV[(1 - ib)zGz + (1 + ib)zGz]}dtdx 
<5 [ C'^(pe^\Vz\'^dtdx 

JQo 



+ - 



c 



/ 9^\Gz\''dtdx + \^fx'' [ ^W\z\'dtdx 

JQo JQo 



Now, we choose ^ = 1. By (4.22), we conclude that 

/ ^e^lVzl^dtdx <C(l + b^)\-^ f e^\Gz\^dtdx + X^iJi^ I ^^9 
JQo ^ya^ J On J On 



2| ^|2 



z\ dtdx 
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Finally, combining (4.19) and (4.23), we arrive 

r r (4-24) 

< C(l + b'^)\ e^\Gz\^dtdx + \^n^ / i/e'^\z\^dtdx] , 

'-JQ J{Q,T)xoj 

which gives the proof of Theorem 4.1. □ 

5 Proof of the main results 

In this section, we will give the proofs of Theorem 1.1 and 1.2. Thanks to the classical dual 
argument and the fixed point technique, proceeding as [4], Theorem 1.2 is a consequence of 
the observability result (1.7). Therefore, we only give here a brief proof of Theorem 1.1. 

Proof of Theorem 1.1. We apply Theorem 4.1 to system (1.5). RecaUing that q{-) e 
L°°(0,T;L"(Q)) and using (4.2), we obtain 



AV / ^'^e'^\z\^dtdx + \ix^ I ipe^lVzl'^dtdx 

JQ JQ 

< C(l + b^) \ [ e^\qz\^dtdx + AV^ / (p^e^\z\^dtdx] (5 1) 

'-JQ J{0,T)xuj -I ^ ' 

< C(l + b^){l + \r\^) \ f e^\z\^dtdx + aV / ^^e^\z\^dtdx . 

'-JQ J(0,T)XLU J 



IQ J{0,T)xui 

Choosing A > C(l + 6^)(1 + |r|^) and large enough, from (5.1), one deduces that 

/ <f^e^\z\^dtdx <c [ ^^e^\z\^dtdx. (5.2) 

JQ J{0,T)xuj 

Finally, by (5.2) and applying the usual energy estimate to system (1.5), we conclude that 
inequahty (1.7) holds, with the observability constant C given by (1.8), which completes the 
proof of Theorem 1.1. □ 
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